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Abstract. We propose a new approximation scheme within equation of motion 
approach (EOM) to spin polarized transport through a quantum dot coupled to 
ferromagnetic leads. It has some advantages over a widely used in the literature 
standard EOM technique, in particular when we are interested in spin polarized 
quantities. Namely, it gives the values of the dot spin polarization which are closer to 
the ones obtained within numerical renormalization group (NRG), than the standard 
EOM approach. While restoring the Kondo effect, the spin polarization vanishes 
and the transport becomes unpolarized, in agreement with NRG and a real time 
diagrammatic calculations. The standard EOM procedure gives nonzero values of the 
spin polarization, and the transport is still spin polarized. Both approximations give 
the same correct splitting of the Kondo peaks due to ferromagnetism in the electrodes. 
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1. Introduction 

The growing experimental interest in magnetic materials has evolved into a new research 
field, spintronics, where the transport properties are governed by the electron spin 
rather than charge. Thus, it opened new possibilities for technological progress in 
nanoelectronics [H E] and quantum computing [21 |3]. On the other hand, due to 
continuing experimental progress in miniaturization of electronic devices it became 
possible to study the fundamental problems of quantum mechanics. One of such 
examples is the Kondo effect [1] in quantum dots (QD). 

The Kondo effect can be observed when the dot has unpaired spin or, in other 
words, has odd number of electrons. Thus the unpaired spin on the dot forms a many 
body singlet state with conduction electron spins in the leads. This state manifests itself 
in a resonance at the Fermi energy in the dot density of states and zero-bias maximum in 
differential conductance. The Kondo effect was predicted a long ago [5l[6l[7], extensively 
studied theoretically [8], [9l [THl [H] and confirmed in a series of beautiful experiments 
[T2| [T3l [HI [151 [IS] ill QD coupled to normal (non-magnetic) leads. 

If the normal leads are replaced by ferromagnetic ones, spin degrees of freedom 
start to play significant role modifying transport [T7l[T8l[19|[2nt[2T|[22t[28|[Ml[25l[26l 
[27l|28l[29l[30l|3ll[32l[33l[3l[35]and thermoelectric properties [36], eventually leading 
to new phenomena. One of such new effects, associated with the ferromagnetism in 
the leads, is a sphtting of the Kondo resonance due to spin dependent quantum charge 
fluctuations [20l[26l[27l[28l[33l[3ll[35l[36]. This splitting depends on the magnitude 
of the leads magnetizations as well as on their alignment in both leads. In particular, 
when the leads magnetizations point in opposite directions (anti-parallel alignment), 
one observes no splitting of the resonance and full Kondo effect is present for all spin 
polarizations. However while measuring differential conductance the Kondo resonance 
gets suppressed with increasing of the leads polarization, flnally leading to a complete 
disappearance of zero-bias anomaly for fully polarized electrodes (p=l). In this case 
transport is completely blocked and no current flows for any voltage. On the other 
hand, in parallel alignment, the Kondo resonance is split and also gets suppressed with 
increasing the leads polarization. However, in this case transport is not completely 
blocked even for p=l, as one spin channel is still conducting [37j . 

Another interesting phenomenon is a compensation of the Kondo effect by external 
magnetic fleld {B). As it was mentioned before, at zero magnetic fleld in parallel 
configuration, one observes splitting of the Kondo resonance. Moreover, finite spin 
polarization is induced on the dot due to ferromagnetic electrodes. It turns out, that 
applying external magnetic field one can recover full Kondo effect, i.e. no splitting of 
the zero-energy resonance and vanishing of spin polarization on the dot. The transport 
becomes unpolarized. Thus, at certain magnetic field B = Bcomp, which we call it 
compensating field, the strong coupling limit is recovered. However, the problem is 
with standard equation of motion approach (EOM) [20l [36] , as it gives non-zero spin 
polarization and spin polarized conductance, even at B = Bcomp, i-e. when there is no 
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splitting of the Kondo resonance. This is in contradiction with other approaches, hke 
numerical renormalization group (NRG) [27] and a real-time diagrammatic technique 
[M] , which correctly give zero spin polarization and equal spin dependent contributions 
to the conductance at i? = Bcomp- It is the purpose of the present paper to show how one 
can cure a standard EOM approach from their disabilities of non-zero spin polarization 
and spin polarized conductance at i? = -Bcomp- 

The paper is organized as follows: in Sec. [2]the model and details of calculations are 
presented. Section [3] shows a comparison of density of states and the on-dot occupations, 
obtained in different approaches. Section H] is devoted to compensation of the Kondo 
effect within the standard and the present EOM approaches, and finally conclusions are 
given in Sec. [51 

2. Model and formulation 

Our model system, i.e. quantum dot coupled to external leads is represented by a single 
impurity Anderson model in the limit of strong on-dot Coulomb repulsion {U —>■ oo). 
The model Hamiltonian in the slave boson representation [381 133 ESl HQ] reads 



where X = L (R) denotes left (right) lead, c'l^^^ (cAko-) is the creation (annihilation) 
operator for a conduction electron with the wave vector k, spin a in the lead A, f^ {f„) 
is a fermion operator, creating (annihilating) spin a on the dot, while 6"*" (b) is a boson 
operator responsible for creating (annihilating) an empty dot state. The product of the 
fermion and boson operators gives a real dot electron operator {d^ = fa)- Vw^^ is the 
hybridization between localized electron on the dot with the energy and conduction 
electron of energy eAk in the lead A. Ferromagnetism of the electrodes is modeled via 
spin dependent conduction energy bandwidths. The constraint of no double occupancy 
is exactly taken into account by the non-canonical commutation rules for fermion and 
boson operators [38j . 

The total current / = Jo-, flowing through a quantum dot is given in the 
standard form [H] 



where V\rj{uj) = 27r |VAkcrp5(^ — ^xka) is the couphng parameter, Po-(a;) is the on- 
dot spin dependent density of states, and f\{uj) = f{uj — ^\) is the Fermi distribution 
function in the lead A with the chemical potential ^\ and temperature T. 

In order to get the density of states Paij^) one has to calculate on-dot retarded 
Green function (GF) G^^iu). Within equation of motion approach the resulted GF is 





(1) 



Ak 
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with non-interacting {U = 0) 



Ak 

and interacting self-energy 



= y I^Ak|VA(6Ak-.) ^ ^ 

Ak 

which is responsible for the generation of the Kondo effect. The non-interacting self- 
energy So(t(c(j) is an exact solution of the problem with no Coulomb interactions present. 
The interacting self-energy T,j^{uj) is obtained by neglecting terms in the equation of 
motion for G^^{uj) which cannot be directly projected onto original dot GF at this stage 

In the standard EOM approach the ferromagnetism in the leads is modeled via 
spin dependent coupling parameters Txa{^^) = '^'^l^xlPxaiEp), where Ep is the Fermi 
energy in the lead A. In order to get the splitting of the Kondo resonance due to 
ferromagnetism in the leads, one replaces bare dot energy level 6^ in self-energy S7o-(w) 
(Eq. [5|) by renormalized one e„, self-consistently found from the relation [20] 

e„ = + RelEo^ie^) + J:i^{e^)]. (6) 

In this way obtained splitting of the zero energy resonance remains in good agreement 
with a pour man scaling approach [20j. However such a procedure has very important 
drawback, when the Kondo effect is compensated by external magnetic field B. While 
a.t B = Bcomp there is no splitting of the Kondo resonance, as it follows from other 
approaches [20, EH El], it gives non-zero value of on-dot spin polarization and not equal 
spin polarized contributions to the transport, which is in contradiction with NRG [2Z| 
and the real time diagrammatic calculations j34j . 

To see why the standard EOM approach gives non-zero spin polarization at 
B = Bcomp let us examine the structure of the dot GF (Eq. ([3])). At zero magnetic field 
[si = Si = Eq) the denominator of GF can be written in the form 

u-eo-il+ p)^o{uj) - (1 + p)^i{oo + As) (7) 

for spin up electrons, and 

u-eo-{l- p)^o{uj) - (1 - p)^i{uj - Ae) (8) 

for spin down electrons, respectively. Both T^qIuj) and Tij{uj) contain spin independent 
couplings to the leads, and the spin dependence is shifted to the polarization parameter 
p, and Ae = — is calculated form Eq. IQ. At B = Bcomp those equations are 

u;-e^-il+ p)So(a;) - (1 + p)S,(a;) (9) 



a;-q-(l-p)So(a;)-(l-p)S,(a;) 



(10) 
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where e„ = eQ + uBcomp, and cr ± 1. As one can notice, the real parts of both equations 
give the same values, that means the DOS for both spin directions will have a charge 
fluctuation resonance centered around the same energy. So in fact it should give the 
same occupations n-^ = n^. However, this is not the case, as we have to consider also the 
imaginary part of self-energies ^0(0;) and S/(u;), which both depend on the polarization 
p. This leads to the smaller spin up occupation due to larger broadening of the charge 
fluctuation resonance in the DOS by factor 1 + p. Similarly, is larger due to a factor 
1 — p in the imaginary part of self-energies. As a result, spin polarization is non-zero at 

To cure those inconsistencies of the standard EOM technique, we propose 
modifications of the approach. We start from two requirements: (i) at zero magnetic 
field it should give the same splitting of the Kondo resonance as the standard EOM 
approach, and (ii) at i? = Bcomp it should lead to zero spin polarization and no splitting 
of the Kondo resonance. This could be easily obtained if we skipped l±p factors in the 
imaginary parts of self-energies but left them in real parts. However such a procedure 
seems to be difficult to substantiate. Here, we propose slightly different approach which, 
by the construction, fulfills the above requirements. Namely, we use the same values 
of Tx] = Txi = Txo for both spin directions but different spin dependent bandwidths 
Dxa in the electrodes, which leads to spin asymmetry in the electrodes. This is one of 
the ways of modeling a ferromagnetism in the electrodes, closely related to the Stoner 
model, but not unique [27]. The splitting of the Kondo resonance due to ferromagnetic 
leads {Ae = e-f — e^) is obtained in the same way as in standard EOM approach (Eq. [6]) 
with replaced Txa in and S/^ by Txa = (1 + <^P\)^xo, where px is the polarization 
in the lead A, and cr = ±1. Note that we use the same bare Fao for both spin directions 
in expression for the GF (Eq. ([3])). The polarization p in the electrodes is calculated 
from spin dependent electron concentrations in the leads, i.e. p = — N^. Those 
concentrations will be different due to unequal bandwidths for spin up and spin down 
electrons. In this way obtained splitting of the zero energy resonance is the same as that 
obtained in pour man scaling and standard EOM technique pO]. Thus the requirement 
(i) is fulfilled. 

To fulfill the requirement (ii) let us write down the denominator of the dot GF (Eq. 
([3])), similarly as for standard EOM equations. At zero magnetic field it gives 

w-5o-So(w) -S^(a; + Ae) (11) 
for spin up electrons, and 

c<j-£o-So(^) -S/(^- As) (12) 
for spin down electrons, respectively. At B = Bcomp the above equations read 

cij — £| — T,q{uj) — T,f{uj) (13) 

u-e^-j:o{uj)-J:j{u) (14) 

where again = Eq + aBcomp- As we can see, there is no splitting of the Kondo 
resonance, however, it also leads to non-zero spin polarization due to different real parts 



Compensation of the Kondo effect in quantum dots coupled to ferromagnetic leads ... 6 

of the above equations. In this case the charge fluctuation resonances in spin polarized 
DOS are centered at different energies. Perhaps the easiest way of achieving zero spin 
polarization aX B = Bcomp is to replace the original dot energies = £o + crB in the 
above equations by the same energies Sa- but expressed in terms of e^- and ReT,j{£f^), 
calculated from Eq. IQ. Thus at 5 = we get splitting of the Kondo resonance and 
non-zero spin polarization, while at B = Bcomp the splitting vanishes as well as the spin 
polarization. However, such a procedure is equivalent to leaving the original dot energy 
level £0 unchanged, even in the presence of the external magnetic field. In other words, 
we have to assume that the B field modifies T,j{uj) only - there is no Zeeman splitting 
of the dot energy level. 

Here we propose another, a more natural way of getting the zero spin polarization 
at B = Bcomp- Namely, we replace the dot energies by e^, not only in the interacting 
self-energy S/^r (Eq. [5]), as in standard EOM, but also in full retarded Green function 
G^. (Eq. [3]). This corresponds to the same replacement of the energies in Eqs. (fT3ll 
and dHl), which in turn leads to the same expressions for Green functions for both spin 
directions at B = Bcomp- Such a renormalization of e^, in a heuristic way, represents 
fact that higher order GFs in EOM procedure give also V"^ contributions (similar to 
J]i^{uj)) going like ln{uj) near the Fermi energy [16], which are neglected in standard 
EOM approach. Perhaps this is a simplest way of enhancing spin correlations, i.e. a 
transfer of the spectral weight from charge to spin sector by the renormalization of the 
dot energy level. This in turn leads to better description of the Kondo effect, as it will 
be shown later on (see discussion in Sec. [3]). Such a procedure gives the correct splitting 
of the Kondo resonance Ae, equal spin dependent contributions to the conductance as 
well as the vanishing of the spin polarization at B = Bcomp- Thus both requirements 
(i) and (ii) are automatically fulfilled. Of course, such a modification does not cure all 
the shortcomings of standard EOM approach, like basic Fermi liquid relations, which 
are little bit less but still violated. However, it leads to better qualitative description of 
the Kondo effect in quantum dots. 

In the following we show how those modifications of the approach influence the 
properties of the quantum dot in the presence of ferromagnetism. In numerical 
calculations all the energies are measured with respect to the Fermi energy Ep = 
in units of F = Tlq + Tjiq = 1. 

3. Density of states 

The most influenced quantity by the above modifications of the EOM approach is the 
density of states (DOS) A^{lu) = -^Gl{uj). 

Figure [1] shows a comparison of spin dependent densities of states obtained in 
the standard EOM (top panel) and in the present, modified EOM approach (bottom 
panel). First of all, one can see that the positions of the Kondo resonances remain 
unchanged in both approaches, giving the same value of the splitting. However, there is 
a change in their spectral weights. While in the standard EOM the Kondo resonance for 
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Figure 1. Comparison of spin up (solid line) and spin down (dashed line) dot density of 
states obtained in the standard EOM (top panel) and in the present approach (bottom 
panel). The model parameters are: T=0.005, £cr=-2 in units of F. The polarization of 
the electrodes is pz,=p_R=0.2. 

spin down electrons is higher, in the modified equation of motion technique (MEOM) 
the resonance for spin up electrons is more pronounced, similarly as in real time 
diagrammatic calculations [31]. Moreover, the DOS around the dot energy level e^r 
also changes leading to higher values for spin up electrons. This is also in agreement 
with real time diagrammatic approach (see Fig. 6 of Ref. [31]). At the same time the 
standard EOM gives comparable or only slightly higher DOS around e^. At the first 
sight, the higher Kondo resonance for spin up electrons seems to be counterintuitive, as 
there is more spin up electrons in the leads, thus they should lead to better screening 
of spin down electrons on the dot. However, one has to remember that due to spin 
dependent renormalization of the dot energy level, the spin up electron occupation on 
the dot is larger, thus the cotunneling processes, including Kondo ones, are more efficient 
in spin up channel. 

As one can notice, the MEOM approach seems to give to large renormalization of 
the dot energy level e^, shifting the charge fiuctuation resonance in the DOS towards 
the Fermi energy for both spin directions (compare both panels of Fig. [1]). However 
such a renormalization infiuences also the DOS around Ep leading to better description 
of the Kondo resonance. This can be seen in Fig. [2], where the densities of states for 
unpolarized leads (p=0) obtained in different approaches are shown. The solid line is 
obtained within non-crossing approximation (NCA) [18] , which is a widely accepted and 
reliable technique for description of the Kondo problem in the case of non-magnetic leads 
and in the lack of external magnetic field [IDl HSl EDI [H] , the dashed one within standard 
EOM, while the dotted one within the present approach. Clearly, the present approach 
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Figure 2. The dot density of states obtained in NCA (solid line), standard EOM 
(dashed line) and modified EOM (dotted line) approach. The model parameters are 
the same as in Fig. [l]but the electrodes are unpolarized now (p=0). 



gives better behavior of the DOS in the low energy regime, thus better description of 
the Kondo effect. 

Physically, as it was previously mentioned, the renormalization of in GF 
represents V"^ contributions to GF, which are obtained while calculating higher order 
GFs in EOM procedure. In particular, performing an effective V'^ order EOM 
calculations, one gets V'^ contributions, going like ln{u)) around the Fermi energy 
[16] . Those contributions are very important for the Kondo effect, as they represent 
inelastic scattering processes leading to the broadening of the Kondo resonance (compare 
Fig. [2] and Fig. 1 of Ref. ^6]). Here similar effect is achieved simply just by 
renormalization of the dot energy level, which in turn leads to a transfer of the 
spectral weight form charge to spin sector enhancing the Kondo effect. Furthermore, 
this can be qualitatively explained even within the EOM approach. Namely, while 
calculating original QD Green function one obtains, usually neglected, an energy 
independent term A = — ^^^^ ^k(/-o-f'CAk-CT), which shifts the QD bare energy level 
towards Ep. In equilibrium, the A can be expressed in terms of the dot GF, i.e. 



QD energy level in the same way as we propose here. This contribution has to be also 
calculated self-consistently, as it depends on the QD energy level, and explains why 
we replaced original QD energy level 6^ by Sa also in full retarded QD Green function 



Another important consequence associated with the above modifications of the 
EOM approach is the average dot occupation n = + n^. For unpolarized leads and 
the dot energy level e„ = —2, the MEOM gives umeom = 0.87 which is close to value 
obtained within NCA {rij^cA = 0.85), while the standard EOM deviates from ujqcA by 
almost 10%, giving ueom = 0.93. 




Clearly, this is contribution renormalizing 



Gliiv) (Eq. (ED). 
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4. Recovery of the Kondo effect 

Now, let us discuss the effect of the external magnetic field on the Kondo effect. In 
the following we assume that the external magnetic field acts on the dot spin only and 
disregard its influence on the properties of the leads. In real experiment, this cannot be 
neglected, as it can lead to the modifications of the magnetic properties and the density 
of states in the electrodes [51]. However, it needs fully self-consistent calculations, which 
are out of the scope of the present work. 

It turns out that applying magnetic field one can recover the Kondo effect. At 
certain field B = Bcomp there is no splitting of the zero energy resonances (in real 
experiments the B field will also modify DOS in the electrodes [51]), the transport 
is unpolarized, and the spin polarization should vanish in this case, so the strong 
coupling limit is reached. Thus, full Kondo effect is recovered. Those conclusions 
have been obtained within NRG [27| and a real time diagrammatic technique [31]. Also 
standard EOM approach gives no splitting of the Kondo resonance in this case [20] but 
does not fulfill the conditions of vanishing spin polarization and equal spin dependent 
contributions to the conductance at Bcomp- The proposed modifications of the standard 
EOM substantially improve the results leading also to unpolarized transport and to zero 
spin polarization at Bcomp- 

Figure [3] shows a comparison of the spin polarization — ni vs. magnetic field 
B, calculated within the present approach (solid line) and the standard EOM technique 
(dashed line). The dotted line represents the splitting of the zero energy resonance Ae^r- 

0.50 

0.25 

^ 0.00 
c 

-0.25 
-0.50 

0.0 0.1 0.2 0.3 0.4 

B 

Figure 3. The spin polarization n| — calculated within the present approach (solid 
line) and within standard EOM (dashed line). Note that within MEOM approach 
spin polarization vanishes &i B = Bcomp = 0.11, while EOM gives non-zero value of 
it. The splitting of the zero energy resonance is also shown (dotted line). The model 
parameters are the same as in Fig. [TJ 

It is clearly seen that the spin polarization obtained within MEOM vanishes, while 
standard EOM gives non-zero value of it at S = Bcomp = 0.11, for which A^o- = 0. 
At this field standard EOM gives n-^ — = —0.135. The spin polarization vanishes at 
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much smaller B field, i.e. aX B = 0.042. 

The non-zero spin polarization results form different densities of states for spin up 
and spin down electrons, which can be seen in Fig. H] (dashed and dotted line) or Fig. 
1(d) of Ref. [20]. Clearly, the spin down DOS is much larger than spin up one, even if 

0.35 
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0.00 

-4 -3-2-1012 
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Figure 4. Spin dependent density of states ai B — Bcomp- Modified EOM gives 
exactly the same DOS for both spin directions (solid line), while the standard EOM 
gives different densities of states for spin up (dashed line) and spin down electrons 
(dotted line). The parameters are the same as previously used. 

there is no splitting of the Kondo resonance. Similarly, other approaches, like NRG or 
real time diagrammatic approach, give also asymmetric density of states but at the same 
time they give zero spin polarization. This is achieved in the following way. Spin down 
electrons have larger DOS at the Fermi energy but smaller around e^; in comparison 
with spin up electrons, so the resulting DOS integrated with the Fermi distribution 
functions give the same occupations, thus no spin polarization [271 . The MEOM 
approach also gives zero spin polarization but the price we have to pay for this are the 
same densities of states for both spin directions (see solid line in Fig. H]). 

It is interesting to see the behavior of the dot spin-dependent occupations as a 
function of the leads polarization in both approaches. Corresponding spin up, spin 
down and the total occupations are shown in Fig. [5] at zero magnetic filed (left panels) 
and at B = 0.05 (right panels). The top panels represent occupations obtained with 
help of standard EOM, while bottom panels - those obtained by the present, modified 
EOM approach. As one can see, a.t B = both approaches give n-^ = for unpolarized 
leads (p=0). However, the difference between and n|, (spin polarization) obtained in 
the standard EOM is much smaller than in the modified EOM approach. Unfortunately, 
both approaches give too small values of spin polarizations in comparison to the other 
approaches ^71 |3l] . At finite magnetic field, modified EOM gives = for a finite 
p, for which Ae^j = 0. On the other hand, the standard EOM gives ^ n^, in 
contradiction to NRG results (compare Fig. 1 of Ref. [27]). 

While the standard EOM approach gives worse (smaller) values of spin polarization, 
it gives better behavior of the total occupation {n-^ + n^) vs. leads polarization. In the 
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Figure 5. Spin-dependent occupation of the dot at i? = (left panels) and B ~ Q.l 
(right panels) as a function of spin polarization p. Top panels represent results obtained 
within standard EOM, while bottom ones - within the present approach. At B — 
both approaches give n-\ = ni for p=0. At B — 0.05 MEOM gives n| — for a finite 
p, for which Ae^r — 0, while the EOM gives n-f ^ in this case. 



standard EOM approach the total occupation weakly increases with decreasing of the 
leads polarization p (top panels of Fig. [5]), surprisingly giving a correct position of the 
maximum of the total occupation, in agreement with NRG results [27]. In the present 
approach, the situation is opposite, namely, the total occupation weakly decreases with 
decreasing of the leads polarization, giving a correct position of a minimum, not a 
maximum of the total occupation. In both cases, however, the changes of the total 
occupations due to the leads polarizations are very small. 

Finally, let us turn our attention to the transport properties. Figure [6] shows a 
linear conductance Gun = -j^^lev^o vs. external magnetic field B, calculated within 
the present approach (top panel) and the standard EOM technique (bottom panel). 
While both approaches lead to similar behavior of the total conductance Gun, i.e. a 
maximum at B = B^omp = 0.11 (indicated by thin vertical line), they give different 
behavior of spin polarized contributions G^ to it. Within the MEOM approach (top 
panel) the transport is governed by majority spin electrons at small B fields {B < Bcomp), 
and by minority spin electrons at higher fields [B > Bcomp)- Such asymmetry of G„ 
can be explained by the spin state of QD, similarly as the asymmetry of the DOS (see 
Fig. [1]). The peak in G| (dashed line) steams from the fact that at i? < Bcomp the QD 
is occupied by spin up electrons {n^ > ui) (see Fig. [3]), thus the spin up component of 
the cotunneling current is dominant. At S > B^omp the situation is opposite, namely, 
?7,| < Tij, thus spin down cotunneling current is larger. At B = B^omp both spin channels 
equally contribute to the transport, again, indicating that full Kondo effect is recovered 
in this case. Note that, the calculated spin polarization vanishes in this case (see Fig. |3]). 
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Figure 6. The linear conductance Gun (solid lines) and spin polarized contributions 
to it (dashed and dotted lines) vs. magnetic field calculated within the present 
approach (top panel) and the standard EOM technique (bottom panel). Note different 
behavior of Ga in both approaches. 



Such a behavior of the conductance and the spin polarization remains in good agreement 
with a real time diagrammatic technique (compare Fig. 6 (a) of Ref. [31]). Within the 
MEOM approach, the condition = aX B = B^omp steams from the fact that the 
density of states for both spin directions are the same (see Fig. Hj) and the couplings to 
the leads are equal (F^j = F^j = Fao)- 

The situation is somewhat worrying in the standard EOM approach (bottom panel), 
as it gives an opposite behavior of G^- At i? < Bcomp, larger contribution comes form 
minority electrons, and at i? > Bcomp, it comes from majority electrons. Moreover, at 
B = Bcomp, the transport still is spin polarized (G^ 7^ Gj). Accidentally, the transport 
becomes unpolarized at slightly larger field than Bcomp, while spin polarization vanishes 
at different filed B < Bcomp (compare Fig. [3]). 

Under non-equilibrium conditions, the discrepancies between standard EOM and 
MEOM approaches look similar. Figure [7] shows a comparison of the differential 
conductance G{eV) = ^^Jyy vs. bias voltage eV = — /iR at zero magnetic field 
(top panel) and at B = Bcomp = 0.11 (bottom panel), obtained within the present 
approach (solid lines) and the standard EOM scheme (dashed lines). Both approaches 
give qualitatively similar behavior of the total conductance, i.e. the same splitting 
of the zero bias anomaly at zero magnetic field and its lack when the Kondo effect 
is compensated {B = 0.11). However, the present approach gives larger values of 
the conductance. Moreover, within the present approach (not shown here) both spin 
channels equally contribute to the transport at any voltage, when the Kondo effect is 
recovered by the external magnetic field B = Bcomp- On the other hand, standard EOM 
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Figure 7. Differential conductance G{eV) = ^j^p^ vs. bias voltage eV = /iL — hr 
of a quantum dot coupled to ferromagnetic leads with p=0.2. Top panel represents 
G{eV) without external magnetic field, while the bottom one is for B = Bcomp — 0.11. 
Solid lines are obtained with help of the present (MEOM) approach, and dashed lines 
- within standard EOM technique. 



technique leads to a spin polarized G{eV) (see Fig. 2 (d) of Ref. [20]). 

Corresponding voltage dependence of the spin polarization is shown in Fig. [81 
Again, at zero magnetic field (top panel), both approaches give qualitatively similar 
behavior of the dot spin polarization for not too large voltages, in agreement with a real 
time diagrammatic calculations plj. At higher voltages (not shown), in the standard 
EOM approach the spin polarization weakly increases with eV, while it still decreases 
in the present approach. On the other hand, when the Kondo effect is compensated by 
the external B field (bottom panel), standard EOM approach gives non-zero values of 
the spin polarization, while the modified one gives vanishing of it at eV = and very 
weak increasing with the voltage. 

All the above results indicate that the present approach should better describe the 
compensation of the Kondo effect, as the transport (linear and non-linear conductance) is 
unpolarized, there is no splitting of the Kondo resonance and spin polarization vanishes 
in this case, similarly as in other approaches [271 On the other hand, standard 
EOM scheme also gives no splitting of the Kondo resonance but at the same time it 
gives non-zero spin polarization and spin polarized transport properties. 

5. Conclusions 

In conclusion we have proposed simple modifications of the standard equation of 
motion approach to get a better description of the Kondo effect in a quantum dot 
coupled to ferromagnetic leads. Special emphasis was put to the compensation of the 
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Figure 8. Spin polarization n-^ ~ "-i as a function of applied bias voltage eV in 
zero magnetic field (top panel) and at i? = Bcomp = 0.11 (bottom panel). Solid 
lines represent the results obtained by MEOM, and dashed lines - by standard EOM 
approach. Note that within the standard EOM spin polarization does not vanish at 



Kondo effect by external magnetic field. In particular, the present approach correctly 
gives both, no splitting of the Kondo resonance and the vanishing of the dot spin 
polarization at compensating magnetic field, while the standard equation of motion 
approach gives non-zero value of the polarization. Moreover, the transport also remains 
unpolarized within the present approach, while the standard EOM gives not equal spin 
dependent contributions to the linear and non-linear conductance when the Kondo effect 
is recovered. In spite of lack of reliable techniques for studying the non-equilibrium 
transport, the present approach could help us to understand the properties of a quantum 
dot in the presence of ferromagnatism in the electrodes. 

Acknowledgments 

This work has been supported by the Polish Ministry of Education and Science under 
Grant No. N202 081 31/0372. 

References 

[1] G. A. Prinz, Science 282, 1660 (1998). 

[2] Semiconductor Spintronics and Quantum Computation, ed. D. Awschalom et at, Springer, New 
York (2002). 

[3] D. Loss, D. P. DiVincenzo, Phys. Rev. A57, 120 (1998). 

[4] A. C. Hewson, The Kondo Problem to Heavy Fermions, Cambridge University Press, Cambridge 
(1993). 



Compensation of the Kondo effect in quantum dots coupled to ferromagnetic leads ...15 



L. I. Glazman, M. E. Raikh, JETP Lett. 47, 452 (1988). 
T. K. Ng, P. A. Lee, Phys. Rev. Lett. 61, 1768 (1988). 
A. Kawabata, J. Phys. Soc. Japan 60, 3222 (1991). 

Y. Meir, N. S. Wingreen, P. A. Lee, Phys. Rev. Lett. 66, 3048 (1991); Phys. Rev. Lett. 70, 2601 
(1993). 

S. Herschfield, J. H. Davies, J. W. Wilkins, Phys. Rev. Lett. 67, 3720 (1991); Phys. Rev. B46, 
7046 (1992). 

N. S. Wingreen, Y. Meir, Phys. Rev. B49, 11 040 (1994). 

M. Krawiec, K. L Wysokinski, Phys. Rev. B66, 165408 (2002). 

D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. Abusch-Magder, U. Meirav, M. A. Kastner, 

Nature 391, 156 (1998). 
S. M. Cronenwett, T. H. Oosterkamp, L. P. Kouwenhoven, Science 281, 540 (1998). 
J. Schmid, J. Weis, K. Eberl, K. von Khtzing, Physica B 256-258, 182 (1998); Phys. Rev. Lett. 

84, 5824 (2000). 

F. Simmel, R. H. Bhck, J. P. Kotthaus, W. Wegscheider, M. Bichler, Phys. Rev. Lett. 83, 804 
(1999). 

S. Sasaki, S. De Franceschi, J. M. Elzerman, W. G. van der Wiel, M. Eto, S. Tarucha, L. P. 

Kouwenhoven, Nature 405, 764 (2000). 
A. N. Pasupathy, R. C. Bialczak, J. Martinek, J. E. Grose, L. A. K. Donev, P. L. McEuen, D. C. 

Ralph, Science 306, 85 (2004). 
J. Nygard, W. F. Koehl, N. Mason, L. DiCarlo, C. M. Marcus, |cond-mat /04104671 



N. Sergueev, Q. F. Sun, H. Guo, B. G. Wang, J. Wang, Phys. Rev. B65, 165303 (2002). 

J. Martinek, Y. Utsumi, H. Imaniura, J. Barnas, S. Maekawa, J. Konig, G. Schon, Phys. Rev. 

Lett. 91, 127203 (2003). 
R. Lu, Z. -R. Liu, |co nd-inat/0210350 [ 

P. Zhang, Q. K. Xue, Y. P. Wang, X. C. Xie, Phys. Rev. Lett. 89, 286803 (2002). 
J. Ma, B. Dong, X. L. Lei, Commun. Theor. Phys. 43, 341 (2005). 
B. R. Bulka, S. Lipihski, Phys. Rev. B67, 024404 (2003). 
R. Lopez, D. Sanchez, Phys. Rev. Lett. 90, 116602 (2003). 

B. Dong, H. L. Cui, S. Y. Liu, X. L. Lei, J. Phys.: Condens. Matt. 15, 8435 (2003). 

J. Martinek, M. Sindel, L. Borda, J. Barnas, J. Konig, G. Schon, J. von Delft, Phys. Rev. Lett. 

91, 247202 (2003). 
M. -S. Choi, D. Sanchez, R. Lopez, Phys. Rev. Lett. 92, 056601 (2004). 
J. Ma, X. L. Lei, Europhys. Lett. 67, 432 (2004). 

J. Konig, J. Martinek, J. Barnas, G. Schon, CFN Lectures on Functional Nano structures , Eds. K. 

Busch et ai, Lecture Notes in Physics 658, Springer, 145 (2005). 
J. Martinek, M. Sindel, L. Borda, J. Barnas, R. Bulla, J. Konig, G. Schon, S. Maekawa, J. von 

Delft, Phys. Rev. B72, 121302(R) (2005). 
Y. Tanaka, N. Kawakami, J. Phys. Soc. Japan 73, 2795 (2004). 
D. Sanchez, R. Lopez, M. -S. Choi, J. Supcrcond. 18, 251 (2005). 

Y. Utsumi, J. Martinek, G. Schon, H. Imaniura, S. Maekawa, Phys. Rev. B71, 245116 (2005). 

R. Swirkowicz, M. Wilczyhski, J. Barnas, J. Phys.: Condens. Matt. 18, 2291 (2006). 

M. Krawiec, K. I. Wysokinski, Phys. Rev. B73, 075307 (2006); Physica B378-380, 933 (2006). 

M. Krawiec, J. Phys.: Condens. Matter 18, 6923 (2006). 

J. C. Le Guillou, E. Ragoucy, Phys. Rev. B52, 2403 (1995). 

M. Krawiec, K. I. Wysokinski, Solid State Commun. 115, 141 (2000). 

M. Krawiec, K. I. Wysokinski, Supercond. Sci. Technol. 17, 103 (2004). 

H. Haug, A. P. Yauho, Quantum Kinetics in Transport and Optics of Semiconductors, Springer, 

Berlin (1996). 
H. Suhl, D. R. Fredkin, Phys. Rev. 131, 1063 (1963). 
Y. Nagaoka, Phys. Rev. 138, All 12 (1965). 



Compensation of the Kondo effect in quantum dots coupled to ferromagnetic leads ...16 



[44] J. A. Appelbaum, D. R. Pcnn, Phys. Rev. 188, 874 (1969). 

[45] C. Lacroix, J. Phys. Fll, 2389 (1981). 

[46] R. C. Monreal, F. Flores, Phys. Rev. B72, 195105 (2005). 

[47] We neglected GFs {{c+^,_^cxi.ah+f-Af^b))^, {{cyk'-aCx).aftMb))^, 
{{Cy]^'a(^\\^-<J.ft„b\f+b))^ and set {{c1^,^,_„Cxk^^b+ U\f+b))^ = ^AA'^kk'/(eAk-a)G^(w), 
in the second iteration of equation of motion for G^{oj) 

[48] H. Keiter, J. C. Kimball, J. Appl. Phys. 42, 1460 (1971); N. Grewe, H. Keiter, Phys. Rev. B24, 
4420 (1981); Y. Kuramoto, Z. Phys. B53, 37 (1983). 

[49] T. A. Costi, J. Kroha, P. Wolflc, Phys. Rev. B53, 1850 (1996). 

[50] M. H. Hettler, J. Kroha, S. Hershfield, Phys. Rev. B58, 5649 (1998). 

[51] M. Pustilnik, Y. Avishai, K. Kikoin, Phys. Rev. Lett. 84, 1756 (2000); M. Eto, Yu. V. Nazarov, 
Phys. Rev. B64, 085322 (2001). 



